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I. INTRODUCTION 



Motivated by the Callan-Rubakov effect in the context of magnetic monopoles [0], studies 
have been carried out recently on the possibihty that cosmic strings can also catalyze baryon- 
number violation with strongly enhanced cross sections. It has been shown that the wave 
function of a fermion scattering off a cosmic string can acquire a large amplification factor 
near the core of the string, leading to enhancement of the processes that violate baryon 
number inside the string [0,0. The catalysis processes that have been studied include those 
mediated by scalar fields and by the grand-unified X and Y gauge bosons in the string core. 
Although strings, in contrast to monopoles, have no magnetic fields outside, fermions can 
interact quantum-mechanically with the long-range gauge fields via the Aharonov-Bohm 
effect. Depending on the flux of the string and the core model used, the enhanced catalysis 
cross sections (per length) can be of the scale of strong interactions in comparison to the much 
smaller geometrical cross section ~ Aj^^^^, where Aqut ~ 10^^ GeV. In the early universe 
when the density of cosmic strings is high, such processes can play important roles, washing 
out any primordially-generated baryon asymmetry 0], or conceivably even generating the 
baryon to entropy ratio observed today. 

Cosmic strings can be produced during certain phase transitions when a gauge group G 
is broken down to a subgroup H by the vacuum expectation value of some scalar field 0. The 
topological criterion for the existence of a string is a nontrivial fundamental homotopy group 
of the vacuum manifold G/H, denoted by 7ri(G/II). For a connected and simply-connected 
G, the general construction of the scalar field at large distances from the string is given by 



Here r is some generator of G, 6 is the azimuthal angle measured around the string, and g{0) 
and g{2TT) belong to two disconnected pieces of H. In the papers referenced in the previous 
paragraph, the scalar field responsible for the formation of the string is taken to have the 
simple form 0(6') = e*'^^0o = e*^0o- As a result, a non-Abelian string can be modeled by a 
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U(l) vortex, and the scattering of fermions in the background fields of the string is governed 
by the Abehan Dirac equation. In general however, for a given (pQ, the generator r can 
be chosen such that e*'^^0o "twists" around the string in more complicated fashion than a 
phase e*^ times 0o- This gives rise to dynamically different strings which are intrinsically 
non-Abelian |^. One expects the complexity and rich structure of such strings to lead to 
interesting effects on fermions traveling around them. In particular, we will demonstrate in 
this paper that for certain r's, the twisting of (j){6) can result in mixing of lepton and quark 
fields, providing a mechanism for baryon number violations distinct from the processes in 
Abelian strings studied previously. 

Since no strings are formed in the minimal SU(5) model, we choose the gauge group 
SO (10) in this paper as an example of grand unified theories in investigating the B- 
violating process. We will construct string configurations, solve numerically for the undeter- 
mined functions, and study the baryon catalysis in the SO(IO) theory, although we expect 
such processes to occur in other non-Abelian theories as well. In SO (10), stable strings can 
be formed when Spin(lO) — the simply-connected covering group of SO (10) — is broken 
down to SU(5)x^2 by the vacuum expectation value of a Higgs field in the 126 rep- 
resentation 0]. The generators of SO(IO) transform as the adjoint 45, which transforms 
as 24 + 1 + 10 + 10 under SU(5). The 24 and 1 generate the subgroup SU(5) xU(l), 
where the U(l) includes simultaneous rotations in the 1-2, 3-4, 5-6, 7-8, and 9-10 planes. 
We are interested in the generators outside SU(5) because to have noncontractible loops at 
all, g{6) in Eq. (1) has to be outside the unbroken H for some 6. We will refer to the U(l) 
generator as and to any of the other 20 basis generators outside SU(5) as ti; we name 
the associated strings as string- Xau and string- ri, respectively. As we shall see, the scalar 
field of string- Ti causes mixing of leptons and quarks while string- Tau is effectively Abelian 
and no such mixing occurs. Properties of string- Tau such as the string mass per unit length 
P] and its superconducting capability in terms of fermion zero modes |p have been studied. 
We will compare it with string- ri, which will be the main subject of study of this paper. 

In Sec. II, we give more detailed discussion of the Higgs 126 and the breaking of Spin(lO) 
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to SU(5)x2^2, and elaborate on the B- violating mechanism due to the nontrivial winding of 
the Higgs field. In Sec. Ill, we write down an ansatz for the field configuration of each string 
and derive the corresponding equations of motion. The numerical solutions and the energy 
of the strings are presented in Sec. IV, where we find that ri-strings have lower energy than 
Tall-strings, probably for the entire range of the parameters in the theory. Having shown 
that such strings are energetically favorable, we turn to the scattering problem in Sec. V, 
where the Dirac equation in the background fields of the strings is solved, and the differential 
cross section for the B-violating processes in string- ti is calculated. We also comment on 
the role of the self-adjoint parameters and compute their values using our string solutions. 
To establish a common notation and to facilitate reading of this paper, we include in the 
Appendix a discussion about the relevant aspects of the spinor representation 16 of SO (10), 
which accommodates a single generation of left-handed fermions. 

II. SO(IO) STRINGS 

There is considerable freedom in the breakings of SO(IO) down to the low energy 
gauge group SU(3) xU(l). Two commonly studied examples include the breaking via 
an intermediate SU(5), SO(IO)— s>SU(5), and the one via an intermediate Pati-Salam 
SU(4)xSU(2)j;^xSU(2)^ ||10|. Details of the symmetry breaking patterns and the Higgs 
fields inducing the breakings can be found in Ref. 6 and the papers by Slansky and Rajpoot 
| rT| . Kibble, Lazarides and Shafi argued that the strings formed during the phase tran- 
sition SO(IO) ^SU(4)xSU(2)lxSU(2)/j become boundaries of domain walls 0. Thus in 
this paper we choose the SU(5) breaking pattern instead for its simplicity. More precisely, 
we study strings formed when Spin(10)^SU(5) XZ2 by the vacuum expectation value of a 
Higgs 126 (j). The nontrivial element of Z2 corresponds to rotation by 2ti in SO(IO). The 
homotopy group 7ri(Spin(10)/SU(5) x Z2) is Z2 ; therefore a Z2 string is formed during this 
phase transition. The subsequent symmetry breakings can be implemented by the adjoint 
45 of SO(IO) and the fundamental 10 in the usual fashion: 
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Spin(lO) ^ SU(5) X Z2 

^ SU(3) X SU(2) X U(l) X Z2 

^ SU(3) X U(l)eni x . (2) 

This Zi string survives all the symmetry breakings since Zi is preserved at low energies. 

The 126 representation consists of fifth-rank anti-symmetric tensors satisfying the self- 
duality condition 

The component which acquires an expectation value (0) transforms as an SU(5) singlet, and 
to write it down explicitly, we first specify how the SU(5) subgroup is embedded in SO(IO). 
The fundamental representation of SO(IO) consists of 10x10 matrices, which can be labeled 
by indices i, i = 1, . . . , 10 . The generators of SO(IO) in this representation can be written 
as antisymmetric, purely imaginary matrices. The generators of SU(5) in the fundamental 
representation are hermitian, traceless 5x5 matrices which can be written as 

7"a/3 = So,^ lAo,^ , (4) 

where a,/? = 1,..,5 label the matrix elements, and S^A are real 5x5 matrices, repre- 
senting the real and imaginary parts of r. Hermiticity and tracelessness of r require 
'S'a/3 = 5'/3a,^a/3 = " ^/3a, and TtS = 0. A natural way to embed SU(5) in SO(IO) is 
to treat fivc-dimensional complex vectors as ten-dimensional real vectors, i.e. replace the 
paired indices (a, a), where a — 1, . . . , 5 label a five-dimensional vector and a = 1, 2 label its 
real and imaginary parts, by the index i, i — 1,...,10. Then, the generators of the subgroup 
SU(5) of SO (10) can be expressed as 

i{AapIab + SajsMab) , (5) 

where / is the 2x2 identify matrix and M — ia2 , (J2 being the second 2x2 Pauli matrix. 

One can convince oneself that in this (a, a) notation, the rank-five antisymmetric Levi-Civita 
tensor ea^a2a3a4,az which transforms as an SU(5) singlet in the SU(5) notation becomes 



Q;ia20!30!40!B ' 



(6) 



where /(oi ... 05) is defined to equal the number of that takes the value 2. It is also 
straightforward to check that this expression satisfies the self-duality condition (Eq. (3)). 
Thus {(f)) is written as 



where // is a parameter. 

Some words about our notation. The tensor indices ii, . . . ,{5 of 0ii...ig will be suppressed 

for convenience and legibility whenever no ambiguity should arise. In the expressions like 
t4> and e*'^^0 where r operates on 0, r is understood to be in the same representation of (f), 
i.e. T is the short-hand for 

With the symmetry breaking Spin(10)^SU(5)x22; strings are formed. At spatial in- 
finity, the general form of (f) is given by Eq. (1). For the energy to be finite, the covariant 
derivative of 0, = -|- 6^4^^ , has to vanish at spatial infinity; therefore the gauge 
field takes the form = i^r, A'' = , as r —> oo. In the core of the string, there is a 
magnetic flux § A ■ dl — pointing in the direction of r in group space. Strings carrying 
flux pointing in different directions in group space are topologically equivalent since the only 
nontrivial winding number here is one, but dynamically they can differ. Because the scalar 
field (f){9) varies with 9, the embedding of the unbroken subgroup SU(5) in SO(IO) outside 
the string also varies with 9. More precisely, the generators , a = 1, . . . , 24 of the unbroken 
SU(5) at 9 are related to the generators Tq of the unbroken SU(5) at 6* = by the similarity 
transformation 



Consequently, the fermion fields which transform as 1, 5 and 10 under SU(5) are also rotated 
as one goes around the string. How the fields mix depends on which direction in group space 
(f){9) winds. 




(7) 



(9) 
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The SO(IO) generators can be written as 10x10 matrices of the form {T°'^)ij — —i{6f5j — 
S^Sj) , where a, b label the group indices, i,j label the matrix elements, and a, b, i,j all run 
from 1 to 10. In this notation Taii is given by 

r.u^W' + T'' + ... + T'''), (10) 
o 

where the factor of 1/5 is included for (f){6) to have a 27r rotational period. It takes a little 

more effort to write down the ri's. Let us first write the SU(5) generators specified by 

Eq. (5) in terms of r"^ given above. The four diagonal generators are trivial. For the other 

twenty generators, one can group the 10x10 space into 2x2 blocks, and write the 45 r"^'s 
as -7-2a-i, 2/3-1^ ^2a-i, 2/3^ ^2a, 2/3-1 ^2a,2/3^ where a, (3 both ruu from 1 to 5. Then it is not 

hard to see that the twenty linear combinations 

i^^2a-l,2^ _ ^2a, 2/3-1^ ^ 

1 (^2a-l, 2/3-1 ^^2a, 2/3)^ « < /3 (11) 

are all of the form of Eq. (5), and therefore can be chosen to be the twenty off-diagonal 
generators of SU(5). Note that the superscripts a, /3 above label the group indices while the 
subscripts a, /9 in Eq. (5) label the matrix elements. The twenty ti's outside SU(5) then can 
be expressed by the other twenty linear combinations as 

^^ = l(r2a-l, 2/3 + ^2., 2/3-1)^ 

l(^2a-l,2/3-l_^2a,2^)^ « < /3 . (12) 

Other than the SU(5) group properties, the linear combinations above can also be clas- 
sified under the group SO (4), which is locally isomorphic to SU(2) xSU(2). For a given a 
and /3 where a < the two generators of Eq. (11) plus the diagonal 

1j-^2q-1, 2/3-1 _ ^2q,2/3^ 

can be easily shown to obey the SU(2) algebra. Similarly, the two generators of Eq. (12) 
plus 



^ (^^2a-l, 2/3-1 ^^2a, 2/3^ ^^^^j 

generate another SU(2). Thus, for a given a and /9 (a < the six generators of Eqs. (11-14) 
generate rotations in the 4-dimensional space spanned by vectors in the 2q; — 1,2q;,2/3 — 1,2/? 
directions. 



III. FIELD CONFIGURATIONS 

The relevant part of the Lagrangian for the SO (10) theory is given by 

L = ^trF^.F'^'' + {D^<py{Dy) - V{<P) (15) 

where F^^ = -iF^^r„,^^ = -i^^T„ , F^^ = du,A^ - dyA,, + e[yl^, yl^] , = + eA^ ; 
yl^, a — 1, . . . , 45, are the SO(IO) gauge fields and (f) is the Higgs 126. The most general 
gauge-invariant and renormalizable potential V{(j)) contains all the distinct contractions of 
two and four 0's: 

V{(t)) = 'yi0n...i5C...i5 +'^2(0il...i5C...75)^ 

+ '^30nn2n3n4nB 0jin2n3n4n5^W ^2^3^44 ^ji^2€3^44 
+ '^^40ni2n3n4n5 0jij2n3n4n5^n«24^44^ju'24^44 
+ '^50nj2n3n4n5^jii2n3n4n5^n»24^44'/'jij2^3^44 
+ '^60ni2j3n4n50jij2i3ri4n6^iii2j3^440jij2j3^44 ■ 

(16) 

In writing down the ^3 through vq terms above, one has to consider two things: (1) the 
possible ways to contract the indices, and (2) which 0's are to be complex conjugated. 
One can deal with (1) without the complication of (2) by adopting an equivalent real 252 
representation for because a complex, self-dual 12 6- dimensional tensor can be thought 
of as a real, 252-dimensional tensor by dropping the self-duality condition and taking the 
real parts of the resulting complex, 252-dimensional tensor. One can see there are only 
four distinct terms and they are terms ^3 through vq in Eq. (16) above. Then when is 
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taken to be complex, two out of the four 0's have to be complex conjugated to make the 
potential real. There are three possibilities: 4>4>*4>4>*, 4>*4>4>4>*, (fxfxf)* (f)* , for each of the four 
contractions 0000 when is real. But after the self-duality condition is applied, one can 
show that only one of the three terms is actually independent. 

The Euler-Lagrange equations of motion for and are given by 

(17) 

= ie{(D^0)*(r"0) - (r«0)*(D»} , (18) 
where a is not summed over, and where a basis has been chosen so that Tr{T"'T^) = 

We construct for string-Taii a solution of the following form: 
Ansatz I : 

= /(r)e^^^>'Vo = /(r)e^Vo , 

A'^i^r^n, (19) 
er 

A^ ^0, 

where 0o = (0) as defined in Eq. (7). The boundary conditions on the functions are 

/(o) = , /(r) "--^ ^, , 

^(0) = 0, g{r)'^l; (20) 

y(0) is minimized at / = /x. Inserting this ansatz into the equations of motion and using 
the relations TaiiTaii0o = 0o a-nd (Taii0o)*(Taii0o) = 0000 — 3840 = N , we obtain two coupled 
differential equations for /(r) and g{r): 

f" + -f-^^-^f = fiv, + 2Nv,n, 
Tr{r^n) (/ - ^^') = -^Ne'{l - g)f , (21) 
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where the prime denotes differentiation with respect to r, and Tr{r^ii) — | from Eq. (10). 
An expansion of /(r) and g{r) in powers of r around the origin reveals that /(r) is odd in 
r with a linear leading term, whereas g{r) is even in r with a quadratic leading term. 
Inserting Ansatz I for string- Taii into the Lagrangian gives 

_ /'all _ ^M^^l) /2 I ]\Tf/2 I iKj i^ ~ 9^ r2 

+7V(^if + 7V^2/'). (22) 

As a consistency check, note that the equations of motion obtained by varying C^^ with 
respect to the functions g and / are identical to those in Eq. (21). 

Note that the parameters ^3 through ve in the potential V are absent from Eq. (21) and 
above. This is because whenever one index of a given (p is contracted with one index of 
another (p, this index is summed over from 1 through 10, or in the {a, a) notation discussed 
earlier, from a = 1 through 5 and a = 1,2. For a given a, the term with a = 2 by definition 
has an extra factor oi — —1 compared to the term with a — 1. These two terms cancel 
each other when they are added. Because this is true for every a, the third through the 
sixth terms in V vanish identically for the string-Tau ansatz. 

To construct an ansatz for string- ti, we need to consider separately the two sets of 
generators in Eq. (12), which will be referred to as 

ri_ = - (r^"-!' 2/3-1 _ ^2a, 2/3) ^ ^ ^ _ (23) 

As we shall see, it is sufficient to derive the equations of motion for an ansatz based on a 
generator of the form ti+. By a simple redefinition, it will then be possible to construct an 
ansatz based on a generator of the form ti_. For now, we consider the case when ri has 
the form ti+. The simple extension of Ansatz / with Tau replaced by ri does not work for 
string-Ti. The problem arises from the term TiTi0 on the left-hand side of Eq. (17) in which 
a new tensor 0q , 

TlTi0o = 0o, (24) 
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is generated, where 



4>oii...i5 ) if two indices take the values 
(2q;- 1,2/3- 1) or (2q;, 2/3) , 
, otherwise. 



(25) 



As a result, the differential equations for g{r) and /(r) are satisfied only if g{r) = 1 or 
/(r) =0 everywhere, which is not consistent with the boundary conditions given by Eq. (20). 
(Note that the solution g = I and / = is the vacuum field configuration expressed in a 
singular gauge.) 

We construct a nontrivial solution for string-ri by replacing /(r)0o and Tau in Ansatz I 
with (/i(r)0o + /2(?')0o ) respectively. Note that 4>q is not orthogonal to 4>q because 

^oh... 15^011... is 7^ 0- Therefore instead of expanding (f) in 0o and 0q , we will use the more 
convenient basis 0q and 0^ where 



00^ 



^-0-00 



(26) 



and (pQ is orthogonal to 0q : 



^^l...^5*^0^l...^5 



0. 



(27) 



Prom the definition of 0q (Eq. (25)) and the properties of 0o, one can see that 



•«5 



0on...iB , if two indices take the values 
{2a -1,2(5) or (2a, 2/3- 1), 
, otherwise 



(28) 



and (pQ is annihilated by ri: 



Tl(/.^ = 0. 



(29) 



The solution constructed for string-Ti is 
Ansatz II : 
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^ ^^Mr, , (30) 
er 

where as will become clear in the next two paragraphs, the functions /o(r) and /e(r) are 
named after their odd and even parities in r. 

At the origin, we require the fields to be regular. Since 0^ is left invariant by e*"^^^ 
(Eq. (29)) but 0o is not, at the origin /e(0) can be any constant but /o(0) has to vanish. At 
large r, the scalar field (f) has to take the form 

e^- Vo = + 0^) (31) 

for the unbroken gauge group to be SU(5), so both /o(r) and /e(r) approach at large r. 
The boundary conditions on the functions are 

/„(o) = o, /oW'^/x, 

/e(0) = ao, /e(r)™/x, 

^(0)^0, givY^l, (32) 

where ao is a constant. 

The equations of motion for and are closed when the fields take the form in Ansatz 
II . We obtain three coupled differential equations for fo{r),fe{r) and g{r). The algebra 
involved in extracting these three equations, however, is considerably more tedious than in 
the Tall case mainly because the forms of 0q , (pQ and Ti are less symmetric. We will not 
present the algebra involved and simply quote the results: 

■fo + ~fo fo 



fo [v, + Nv,{fl + fl) + y/^M - fe)] 



Tr{Tl)(g"--g']=-Ne\l-g)fl, (33) 



r 
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where e^As = V3 + ^ + ^ + and Tr(Tf) — 1 from Eq. (12). An expansion of g, fo and /g 
in powers of r around the origin gives 

fo{r) = air + agr^ + . . . , 
fe{r) = oo + aar^ + . . . , 

g{r)^b2r'' + hr^ + ... , (34) 

where the coefficients of all the higher terms are related to qq, ai and 62 recursively. The 
function /„ is indeed odd and /g even in r as claimed earlier. 

Inserting Ansatz II for string-ri into the Lagrangian gives 

~ '2^^ + T + /o J + y ^2 /o + Kn. (35) 

where 

Kn. = yfl(/„^ + /e^) + y.2(/„^ + /eT 

+^/Ufo-fl)^- (36) 

Here again, note that the equations of motion obtained by varying with respect to the 
functions fo and /e are identical to those in Eq. (33). 

Now let us consider the other case when ti has the form of Ti_. One can show that 
Eq. (24) now is tiTi^o = 0f , and instead of Ti^f = 0, one has ti0q = 0. Therefore by 
switching the definitions of (pQ and 0^ in Eqs. (25) and (28), all the equations between (24) 
and (32) are preserved, and one can show that the equations of motion are unchanged. We 
conclude that Ansatz //applies to all twenty Ti's, where for ti_|_, 0o and (f)^ are defined by 
Eqs. (25) and (28) respectively, but for ri_, the definitions of the two are reversed. The 
equations of motion are given by Eq. (33) for all cases. 



IV. NUMERICAL CALCULATIONS 



In this section we present the numerical solutions to the two sets of differential equations 
(21) and (33) with the appropriate boundary conditions at the origin and some large value 
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of r. We implemented two methods: the "shooting" and the relaxation methods to handle 



this two-point boundary value problem. In the "shooting" method |]I2[, an initial guess for 
the free parameters at r = was made and then the equations were integrated out to large r 
where the boundary conditions were specified. As the name of the method suggests, the true 
solutions were found by adjusting the parameters at r = in the beginning of each iteration 
to reduce the discrepancies from the desired boundary conditions at large r computed in 
the previous iteration. For string- ri, the small-r expansion of the functions in Eq. (34) gives 
^(0) = ^'(0) = , /o(0) = /^'(O) = /^(O) = , and f'J{0) = la^ , where as is related to ao, ai 
and 62, but the values of 

/e(0) = ao , 
/:(0) = a^ , 

g"{0) = 2h , (37) 

were adjusted to match the boundary conditions at large r. For string-Tau, we have shown 
that /(r) is odd and g{r) is even in r, with /(r) = ar + . . . and g{r) = br"^ -|- . . .. Thus 
only the two values f'{0),g"{0) were free parameters. At large r, discrepancies from the 
boundary condition were corrected by the multi-dimensional Newton- Raphson method which 
computed the corrections to the initial parameters. With an initial guess for the parameters 
at r = 0, this "shooting" process was iterated until the "targets" were met. The fourth-order 
Runge-Kutta method was used to integrate the equations. 

We have also implemented a relaxation scheme for comparison. In this method the first 
step is to express the string energy as a function of the values of the functions / and g (or 
fe, fo, and g) defined on an evenly spaced mesh of points. While a Simpson's rule approxi- 
mation worked well for the middle range of parameters, a more sophisticated approximation 
was used to extend the range of parameters that could be treated. For each interval of 
two lattice spacings, smooth functions / and g were defined by 2nd order polynomial in- 
terpolation from the three mesh points (midpoint and two end points); with the help of a 
symbolic integration program, the integral defining the energy was carried out exactly for 
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the interpolated functions. (By this method the energy obtained is a rigorous upper hmit 
on the true ground state string energy.) To avoid divergences caused by the exphcit factors 
of in the energy density, the first interval had to be treated more carefully — instead of 
fitting the functions with a 2nd order polynomial, we fitted the coefficients of the analyti- 
cally determined power series, such as Eq. (34). Trial functions / and g were chosen, and 
then the energy was minimized by varying each mesh point one at a time, successively going 
through the lattice many times. We found it efficient to begin with a coarse mesh which was 
made successively finer by factors of 2, interpolating the solution at each stage to obtain the 
first trial solution for the next stage. For the final run in each case we used 2048 points. 

We found the results by the two methods to agree to approximately one part in a million 
or better. In general we were able to explore a wider parameter range with the relaxation 
method than with the "shooting" method, but the qualitative features given by the "shoot- 
ing" method remained the same. (The author wishes to thank Alan Guth for implementing 
the relaxation part of the calculations.) 

The dependence of the equations on the parameters in the theory can be simplified if 
r, /, fo and /g are rescaled as {vi < 0) 



Then only the following combinations of parameters appear in the differential equations: 



The Hamiltonian densities 7Y^" and Ti^ for the two strings are simply —C^^^ and —C^ given 
by Eqs. (22) and (35) because all fields are assumed to be time-independent. With the same 
rescaling, one obtains 



r 




(38) 




(39) 




(40) 



15 



and 




2 2r2 



Po+Pe \\ h 



(41) 



2 ; 200A2 



where the equation depends on A2 only but the ti equation depends on both A2 and A3. 

Typical solutions for the two strings calculated from the "shooting" method are shown in 
Figs. 1 and 2, where A2 = 0.132 and A3 = 10.25. For the same A2 and A3, the solutions given 
by the relaxation method appear indistinguishable visually from those in Figs. 1 and 2. For 
string-Tall, we were able to find solutions in the approximate range 10^^ < A2 < 10 using 
the "shooting" method and 10~^ < A2 < 10^ using the relaxation method. For string-ri, 
we explored the range 5 x 10^^ < A2 < 1 and 0.5 < A3 < 10^. In general, the functions 
converged more slowly near the two ends of each range above, and we did not attempt to find 
solutions beyond these limits. We numerically integrated Ti*^" and Ti^ for the solutions we 
computed, and found string- ri to have the lower energy for all the parameters we explored. 
In Fig. 3, the energy density 27TrT-C of the two solutions shown in Figs. 1 and 2 is plotted, 
and the energy of string- Ti is clearly lower. For comparison, we point out that the energy 
per unit length of string-Tau in the range 0.9 < A2 < 4.0 has been calculated by Aryal and 
Everett Our values in this range of parameters agree with theirs to within 1%. 

One of the most important properties of the two strings we investigate in this paper is 
whether string- ri has lower energy than string- Tau. We just showed that this is true for some 
range of the parameters. To systematically explore a wider parameter range, however, it is 
very laborious and time-consuming to calculate the ri solutions for different A2 and A3 first 
and then compute the corresponding energy. Instead, we employ an upper-bound argument 
to reduce the two-dimensional parameter space (A2, A3) to one. We set = fe = fi in 
the Lagrangian and take g{r), f\{r) as trial functions for string- ti. The advantage in using 
fo = fe is that the last term in Eq. (41) vanishes, and the equations no longer depend on 
A3. Moreover, Eqs. (40) and (41) then have the same functional form, differing only in the 
coefficients of the first and the third terms, and one can solve the equations for string- ri the 



16 



same way as for string- Taii using different values of A2. The corresponding energy, denoted 
by Ei{fo = /e), gives an upper bound on the true energy of string- ri by the variational 
principle. If -Ei(/o = /e) < -E'aii for a given A2, then one can conclude that string- ri has the 
lower energy for that value of A2 and all values of A3. (Note that in the limit of A3 ^ 00, 
the trial functions approach the true string solution because for the energy to be finite, the 
last term in Eq. (41) requires fo — > /e-) Our result is presented in Fig. 4, where the ratio 
Eiifo = fe)/E^\\ is plotted as a function of logA2 for 10""^ < A2 < 2.5 x 10^. Note that 
Ei{fo = fe)/Ea.\\ < 1 for all 7 decades of A2, and is approaching an asymptote of 1 (or 
possibly less than 1) as A2 — > 0. For large A2, we find the individual curves of -Eaii vs. logA2 
and El vs. logA2 approach straight lines, suggesting that the ratio -E'i(/o = /e)/-E'aii levels 
off at a constant for large A2. We conclude that string- ri has lower energy than string- 
for 10^^ < A2 < 2.5 X 10^ and all A3, and probably is the ground state for the entire range 
of the parameters in the theory. 



V. SCATTERING SOLUTIONS 

To study the scattering of fermions by an SO (10) cosmic string, one first needs to under- 
stand the 16-dimensional spinor representation of SO(IO) to which the left-handed fermions 



are assigned. Spinor representations certainly have been discussed in the literature |]13 
but to establish a common notation, we discuss in the Appendix the construction of the 
generators, the sixteen states and the identification of states with fermions that are relevant 
to this paper. 

Now we proceed to study the Dirac equation 

{i ^-e A^r" - m)^ = (42) 

in the background fields of string- Tau and ri: A^r" = Aj^'Vaii and Aj^Ti. As shown in the 
Appendix, the fermion fields can be written as a 16-dimensional column vector where each 
component is identified with a fermion given by Eq. (A. 16). The generators Tau and ti can 
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be written as 16x16 hermitian matrices, where Taii is diagonal with one diagonal entry equal 
to |, ten entries equal to ^ and five entries equal to — ^. For ti, we choose ri = |(t^^ + t^^) 
for illustration. We find that Ti takes the block- diagonal form 

\ 

(43) 



-n 



B 
B 



where 



B 



'^0 7^ 





/ 



(44) 



and / is the 2x2 identity matrix. For string-Tan, since Taii is diagonal, Eq. (42) decouples 
into sixteen equations, one for each component of the wave function, and there is no mixing 
of leptons and quarks due to twisting of the Higgs. However, since the sixteen eigenvalues of 
Tall are all fractional, all sixteen fermions scatter nontrivially off the string via the Aharonov- 
Bohm effect. As pointed out by previous studies, the wave functions of these fermions can 
be strongly enhanced near the core of the string, leading to strong B-violating processes 
inside the string. 

In the case of string-ri, upon diagonalizing ti by a unitary matrix U and simultaneously 
rotating the fermion basis in Eq. (A. 16) to ip = Uip, we can write ip as 

ip ^ {e" + u'{,e~ — ul,u'^ + di — di ,U2,ul,ds,d2, 
U3 + d2,U3 - d2,U2 + dl,U2 - dl,ui,u ,e'^ , dDi 



(45) 



and Eq. (42) again decouples into sixteen equations of the form 



(i ^ + e\i 4^ - m)^i = , 



(46) 
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where each ipi interacts with the gauge field with couphng strength eAj ; Aj are the eigenvalues 
of — Ti . The eigenvalues are Aj = ^ for e + , z/"^ + rfi , M3 + ^2 , M2 + c^s , A j = — | for 
e — ul ,1/'^ — di — d2 ,U2 — dl , and Aj = for all others. Since the e + u'^ and e — u*^ 
components have opposite eigenvalues, we expect a pure e or u'^ to turn into a mixture of e 
and u"^ as it propagates around the string, producing baryon-number violation. 

Before calculating the scattering amplitude, we first comment on the choice of gauge in 
this problem. The fields in Ansatz II (See Eq. (30)) for string-ri were constructed in a gauge 
where the scalar field winds with 9 and the gauge field falls off as at large r. The 
particle content, however, is probably most transparent in a different gauge where 4> does 
not wind with 9 and — > at large r everywhere except on a sheet of singularities at ^ = 0. 
We will refer to the former as the 1/r- gauge and the latter as the "sheet" gauge, in analogy 
with the "string" gauge of a magnetic monopole. Continuing to work in the diagonalized 
basis, the fermion fields in the "sheet" gauge, ■0o, are related to those in the 1/r-gauge, tp, 
by the gauge transformation 



3— m(7r- 



(47) 



We will solve the Dirac equation and calculate the scattering amplitude in the 1/r-gauge, 
and then write down the baryon-number violating cross section in the "sheet" gauge. 

In the presence of an infinitely-thin Ti-string along the z-axis, the gauge field Aj^ takes 
the form A^^ — A^^ = 0,^4^^ = ^ , where {r,9) denote the usual polar coordinates with 
9 running counter-clockwise from the positive x-axis. Owing to the symmetry along the 
2;-axis, the matrix 73 in Eq. (46) drops out, and with the choice for the 7-matrices 



7o 



72 



^3 ^ 



-as 



J 



^ -iai ^ 



V 



71 



73 



i(T2 








-i(J2 





A 


-1 





(48) 



i(7i 

Eq. (46) decouples into two independent equations for the upper and lower 2-component 
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spinors of ipi, where the two equations differ by the sign of the mass term. Writing the 
upper spinor of ipi as 



^ XI (r) ^ 



X2{r)e" 



And—iEt 



(49) 



one can show 



m — E —i (ji)r + 
-i [dr - ^) -m-E 



= 0, 



(50) 



VX2y 



and the solutions are Bessel functions of order (n + Aj) and — (n + A^): 



(51) 



The appropriate boundary conditions to impose, as pointed out in Ref. 14, are the square- 
integrabihty of the wave functions near the origin and a self-adjoint Hamiltonian. The usual 
requirement that wave functions be regular at the origin is sometimes too strong and has 
to be relaxed. Since Ji,{r) ~ r'^/{2''u\) for small r, one can see that the solutions above 
are square-integrable if the + sign is chosen for the modes n + Aj > 0, and the — sign for 
n + Xi < —1. For the mode — 1 < n + Aj < 0, however, both choices are square-integrable 
albeit neither is regular at the origin, and the solution takes the form 
/ . . . . \ 

(52) 



sin /X Jn+\i + cos II J-(n+\i) 



\ :BT^^(sin/i Jn+Ai+1 - COS/i J_(„+A,+l)) j 

where /i is the self-adjoint parameter. 

The scattering amplitude f^^{9) for the ith fermion in ip appears in the asymptotic wave 
function written as the sum of the incident plane wave and the scattered part: 



(53) 



where ue and ve are given by 
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Ue 



( 1 ^ 



\ E+m J 



1 ^ 



ve 



(54) 



\ E+m^ / 



Expanding e^^^ = e^''^'^°^^ and e^''^ in Bessel functions using 



ikr cos 9 



ind 



(55) 



and with 



(56) 



Eq. (53) can be matched to the solutions in Eq. (51) mode by mode at large r. Then the 
scattering amplitude can be calculated: 



-i([A,]+i)e 



m /sin(f-7rA, 



sm 



(57) 



where [Aj] denotes the largest integer less than or equal to A^, and 5 is related to Ai and the 
self-adjoint parameter tan /x by |n 



1 — tan /i Ajvr 

tan = tan 

2 



(58) 



1 + tan fi 

With the gauge transformation Eq. (47), one can easily see that {4'o)i in the "sheet" gauge 
is given by Eq. (53) without the phase Q-^^-'i'^-^) , 

To illustrate the processes that violate the baryon number, we consider an incident beam 
of electrons propagating in the fields of the string. We will study the (e, u'^)-subspace and 
ignore other fermions since e in is mixed with u'^ only. In the "sheet" gauge, the eigenstates 
of Ti can be written as 



(A 



e — u 



(59) 



and the electron is simply given by 



(A 

2 



(60) 
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An incident wave of electrons can be written as 



Oinc 



Ue 



(A 

2 



Jikx-Et) 



(61) 



which scatters into 



sea 





(A 






> 


< fHo) 


2 


+ /-'(«) 


> 








1 


> 



(62) 



Note that the suppressed index on the 2-component spinors ue and ve should not be confused 
with the index associated with the 2-component eigenvectors used here to label the e + u'^ 

and e — u'^ components of the Dirac field. Rewriting iposca above as 



+ 



fm-rm 




i{kr-Et) 



(63) 



one finds that the scattered wave consists of a mixture of electrons and li'^-quarks. 

The differential cross section per unit length for the production of w-quark is defined by 



where J* = ■07* ^0 . Substituting ■0omc and ■0osca into the currents, one obtains 

which can be written out using Eq. (57) as 

1 



(64) 



(65) 



da 



(66) 



The calculation above was done in the limit of zero string width. Now let us examine 
the string core. The structure of the string core is "encoded" in the self-adjoint parameter 
6 (or jj,, related to d by Eq. (58)), which appears in the differential cross section above. 
In general the self-adjoint parameter is determined either from physical properties at the 
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origin or sometimes by symmetry arguments. Since the string solutions have aheady been 
obtained in the previous section, we can find /i by solving Eq. (50) numerically for the mode 
— 1 < n + Xi < 0, using the realistic form g{r)/r for the gauge field computed earlier in 
place of the 1/r in Eq. (50). As we have shown, Aj = ±| for the fermions that scatter 
nontrivially off the Ti-string. Thus the special mode satisfying — l<n + Ai<0 takes the 
value n -\- Xi — — |, where n = —1 for Aj = | and n = for Aj = — |. Recall that in the 
calculation of g{r), the radial distance r was rescaled to the dimensionless \f—v\r {y\ < 0), 
where vi is the quadratic coupling in the Higgs potential in Eq. (16). Rescaling X2 and r by 



E + m 

r yJ-VxT , (67) 



X2 I — ^ — X2, 



and replacing Aj in Eq. (50) by Xig{r), Eq. (50) can be rewritten as 

drXl = ^^^Xl + PX2 

2r 



9.X2 = -^X2-/3xi (68) 



for Xi — l,n — —1, and 



drXi = — ^Xi + PX2 

drX2 = ^^^X2 - Pxi (69) 

for Aj = — |, n = 0. The parameter P is defined by 

f3 = k/V^,, (70) 

and the bars over Xi)X2 are used to distinguish the solutions of A, — — | from those of 
Aj = |. Upon closer inspection of the two sets of equations above, one finds that Eq. (69) 
is in fact identical to Eq. (68) if Xi is identified with X2 and X2 with — Xi- What about the 
boundary conditions at the origin? In Eq. (49), for n = — 1, the upper component depends 
on 6 but the lower component does not, and vice versa for n = 0. Therefore xi and X2 
must vanish at the origin for the solution to be continuous, but X2 and Xi can be nonzero 
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at r = 0. One thus has xi = X2 and X2 — — Xi- Since Eq. (68) is hnear, the value of X2(0) 
can be chosen arbitrarily when integrating the differential equations. 

The self-adjoint parameters /j, for Aj = | and /i for = — | are determined by matching 
the solutions of Eq. (68) to the asymptotic expression in Eq. (52) at some radius r. For 
n + Xi — — |, the Bessel functions in Eq. (52) are simply J_^i, which have the analytic forms 

FY FY 

Ji(x) = \ — sinx, J_i(x) = \ — cosx. (71) 

2 \l TTX 2 V TTX 

Then Eq. (52) leads to the simple expression for n and /2: 

Xi 

Xi 



tan(u + /3r) , 

X2 



tan(/i + /3r) , (72) 
X2 

which can be inverted to give and /i at a given r, using xi and X2 computed from Eq. (68). 
Using Eq. (72) and trigonometric identities, one finds 

+ (73) 

Note that the solutions depend on (5 which appears in Eq. (68), and the quartic couplings 
A2, A3 in the Higgs potential. The parameter j3 defined in Eq. (70) measures the ratio of 
the incident fermion momentum k to the Higgs mass parameter \/—vi, which is of the order 
of GUT energy scale. To put it another way, (3 measures the string width relative to the 
wavelength of the incident fermion. In Fig. 5, we set /3 = 1 and plot /i computed from 
Eq. (72) at a given r for three sets of A2 and A3. The true value of ji is given by the limit 
r — > 00. In Fig. 6, we choose the same set of parameter as in Figs. 1-3: A2 = 0.132 and 
A3 = 10.25; is shown for five values of j3 ranging from 0.1 to 2.0. One can see that as j3 
decreases, i.e. when the wavelength of the fermion becomes large compared to the string 
width, jjt. decreases. 
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VI. CONCLUSIONS 



We constructed two types of strings, string-Tan and string- ti, in the SO(IO) grand unified 
theory. They are topologically equivalent but dynamically difi^erent strings, produced during 
the phase transition Spin(lO) —>■ SU(5) x Z2 in the early universe. String-Tau is efi^ectively 
Abelian, and can catalyze baryon number violation with a strong cross section via grand- 
unified processes inside the string. It has been the subject of study in several recent papers. 
The richer Higgs structure of string- ti, on the other hand, has been shown in this paper 
to induce baryon catalysis by mixing components in the fermion multiplet, turning leptons 
into quarks as they travel around the string. The underlying B-violating mechanism is the 
"twisting" of the scalar field, which leads to different unbroken SU(5) subgroups around the 
string. This mechanism is distinct from the grand-unified processes which can only occur 
inside the string core where the GUT symmetry is restored. 

The corresponding string solutions have been calculated numerically with both the 
"shooting" and the relaxation methods. The energy of both strings was computed. With an 
additional upper bound argument, we found string-ri to have lower energy than string- Tau 
in a wide range of parameters: 10~^ < A2 < 2.5 x 10^ and all A3. The ratio of the upper 
bound on ti energy to the Tau energy increases as A2 decreases, and possibly approaches 
one from below as A2 0. Scattering of fermions in the fields of string-ri has also been 
analyzed, and the B-violating cross section is given by Eq. (66). We conclude that string-ri 
is more stable than string-Taii, and can catalyze baryon decay with strong cross sections via 
the interesting mechanism of Higgs field twisting. 
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APPENDIX 



The generators of S0(2n) in the spinor representation can be constructed from a set of 
2" X 2" hermitian matrices r^"\ a = 1, . . . , 2n , which satisfy the Chfford algebra 



Starting with the two Pauh matrices for n = 1 



ab • 



(A.l) 









r(i) _ 

) 2 — 




—I 




.1 


% 









(A.2) 



one can iteratively build the higher- dimensional r^"+^'' from the T'^^ by 



p{n+l) ^ 



-.(n+l) 
2n+l 



-.(n+l) 
2n+2 



( p(n) Q ^ 

^ 1^ 



^1 0^ 



, a = 1, . . . , 2n 



i 



(A.3) 



One can check that these F matrices satisfy the Clifford algebra. The ^"^^^ generators of 
S0(2n) are constructed by 

1 



= ^[ra,r6] , a,b = l,...,2n 



where Mab satisfy the S0(2n) commutation relations 

[Mab, MJ = -t{SbcMad + - 5acMu - hdMac) • 



(A.4) 



(A.5) 



Thus far, we have used the explicit matrix notation to construct V and M . For conve- 
nience, however, we will use an alternative notation in which each of the 2" x 2" matrices 
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is written as a tensor product of n independent Pauli matrices, each acting on a different 
two-dimensional space. We cfioose tfie convention that the first matrix from the right in the 
tensor product acts on the largest 2x2 block in the matrix notation, while the second from 
the right acts on the next, and so on, with the matrix on the left acting on the smallest 2x2 
block. In this notation, the 10 F's of SO(IO) given by Eq. (A. 3) become 

Fi = (71(73(73(73(73 , r2 = (72(73(73(73(73 , 

^3=1 CnO'sO'sC^S , r4 = / a2a3a3(T3 , 

— I I (71(73(73 , Fq — I I (72(73(73 , 

^ I I I CTi(73 , Tg = / / / (72(73 , 

Tg I I lai, rio = ////(72, (A.6) 

and the 45 generators M can be found accordingly. Furthermore, one can write down the 
five diagonal M's that generate the Cartan sub-algebra: 

Mi2 = ^ a^IIII , 
M34 = ^ la^III , 
M56 = ^ Ila^II , 
Mrs = \ lllo^l , 

Mg^o^^IIIIa^. (A.7) 

The eigenvalues of the five generators above can be used to label the states in the 
spinor representation. Let |ei, . . . , be the eigenvalues of Mi2,...,M9io respectively 
with ej = -|-1 or — 1, and denote the states by 

I 6162636465 ) . (A.8) 

This 32-dimensional representation is reducible to two 16-dimensional irreducible represen- 
tations because there exists a chirality operator 

x = (-0'rir2...rio 

= (73(73(73(73(73 , (A.9) 
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which satisfies the commutation relations 



(A.IO) 



Moreover, 



(A.ll) 



where the eigenvalue Yl^et is +1 or —1 depending on whether the number of spins that are 
down (ej = — 1) is even or odd. 

We assign the sixteen left-handed fermions to the states of positive chirality, i. e. states 
with even number of = —1. The explicit identification of states to fermions can be 
achieved by first breaking the SO (10) 10x10 representation into an upper 6x6 and a lower 
4x4 blocks for the subgroups S0(6) and S0(4), and then embedding SU(3) in S0(6) and 



SU(2) in S0(4). The generators for S0(4) are Mat, a,b = 7, 8, 9, 10, and with the choice 



for the generators of SU(2), one can easily verify that the last two spins in | €162636465 ) label 

the SU(2) states with | H — ), | — h ) labeling the doublets and | + + ), | ) the singlets. 

Similarly, the first three spins in I6162636465) label the SU(3) states with | + ++), | ) 

labeling the singlets, and | + H — ), | — h+) with their permutations labeling the SU(3) 
triplets. One also needs the charge operator Q to make the assignment unique. In SU(5), 
Q = diag {1/3, 1/3, 1/3, 0,-1), which takes the form 




(A.12) 



Q = -(M12 + M34 + Msg) - Mg 10 . 



(A.13) 



In the SO (10) spinor representation. 



Q|ei...65) = {^(61 + 62 + 63) -|} 



61 ... 65 ) . 



(A.14) 



Putting all the above together one obtains 
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+ + + ++) 


= 1/% 


1 + + + 






h ++ ) 


= ui, 


h 




= di 


- H h+ ) 




1 - + - 






H h+) 
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1 H 




= dl 


+-) 


= i^, 




— ^) 


— e" 


+ + -+-) 


= «1 , 


1 + + - 


-+) 


= dl 


H h+-) 


= U2 , 


1 H h 


-+) 


^d2 


- + + +-) 


= U3, 


1 - + + 


-+) 


= ds 



(A. 15) 



Since we already know how to express the generators Mab as matrices, we can write the 
states as a single 32-dimensional column vector which is projected into two 16-dimensional 
vectors of positive and negative chirahty by the operator P± = |(1 ± x)- We find 

■0 = {ly'^ ul U2 u%, dz d2 dl e~ uz U2 ui v dl d^ dl)L- (A. 16) 

In this paper, we studied two types of strings: string-Tan, where Taii is given by Eq. (10), 
and string-Ti, where ti can be any of the generators in Eq. (12). It is easy to see that in 
terms of Mab, '^aii is written as 

Tall = I (Mi2 + M34 + M56 + M78 + Mg 10) , (A. 17) 

o 

and I ei . . . £5 ) is an eigenstate of Tau with eigenvalue ^J2i^i ■ ^or the left-handed fermions 
above, — \ for i/'^, ^ for e"*", u, d, u'^, and — ^ for v, e~, d^. 

To study how Ti act on the fermions, we write Ti+ and Ti_ defined in Eq. (23) as a 
product of five Pauli matrices using Eqs. (A. 4) and (A. 6), and then replace the matrices a\ 
and (72 by the usual raising and lowering operators (T± = |((7i ± 1(72)- One obtains 

^^^^ l(^2a-l,2/3^^2a,2/J-l) 

— I . . . la+as . . . a^a+I . . .1 

+1 ...la^as-.-a^a^I ...I (A.18) 
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and 




^ _('^2a-l,2/3-l _ ^2a,2l3\ 

— I . . . la^as . . . (J^(J_I . . . I 



—I . . . la-as . . . a^cr+I . . . I 



(A.19) 



where a, (3 — l,...5,a < (3, and the two a± matrices in each term occupy the ath and 
/3th positions from the left. Now one can read off from the hst of fermions above which 
particles are mixed by a given ti. For generators of the form ti+, one immediately finds 
that except for the case a — A, (3 — 5, all mix leptons with quarks; when a — 4, (3 — 5, 
the generator mixes (e"^, v'^), {ul, df), (mj, rfg), and (wg, rfg). For generators of the form ti_, 
leptons are mixed with quarks when a = 1, 2, or 3 and (3 = 4 or 5. 
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FIGURES 

FIG. 1. The solution of string-ri, g{r), foir), feif), as a function of dimensionless r for the case 
A2 = 0.132, A3 = 10.25. The function g{r) represents the spatial dependence of the gauge field, and 
fo{f)-,fe{r) represent that of the Higgs field. 

FIG. 2. The solution of string- Tq^/, /(r), as a function of dimensionless r for the same case 
as in Fig. |^. Here g{r) represents the spatial dependence of the gauge field and /(r) that of the 
Higgs field. 

FIG. 3. The radial energy density 2TTr7i{r) (in units of V1/V2) of string-ri and Tau, computed 
from the solutions in Figs. 1 and 2. 

FIG. 4. The ratio of the upper bound on n energy, i?i(/o = /e), over the Taii energy, E'aii, as a 
function of A2. Ei{fo = fe) is calculated by setting fo = fe in the Lagrangian. 

FIG. 5. The self-adjoint parameter ^ computed from Eq. (72) at a given r, for three sets of 
(A2, A3): (0.132, 10.25), (0.264, 20.50) and (0.528, 41.0), where P = 1. The true value of ^ is given 
in the limit r — > 00. 

FIG. 6. The self-adjoint parameter fi computed from Eq. (72) at a given r for different ratios 
of P, where A2 = 0.132, A3 = 10.25. 
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